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Abstract 

We studied the oxygen vacancies(Vb) in mtile Ti02 by using GqWq approximation on top of 
GGA+C/ as a method of choice to improve the gap. Since there is no extensive agreement regarding 
the characteristic of electron localization for Ti02, we examine combined GqWo@GGA+U scheme 
in which both are conceptually one step toward better enumeration of the non locality of exchange- 
correlation potential . Our GqWq@GGA+U results realize and confirm the weak nature of electron 
correlation in rutile Ti02 and shows that the [/-dependence of the energy gap in perfect bulk 
is slightly stronger than in defected sample. In addition, we studied the U- dependency of Vq 
defect states and found that different charged vacancies shows different [/-dependence . While 
the application of GqWq correction would improve the quasiparticle gap and formation energies, 
however the Vq states , in contrast to experiment, remains entangled with the conduction band. 
Finally, we used PBEO and HSE06 hybrid functionals and found that these exchange-correlation 
functional particularly HSE06 that reproduced the real gap and provide desire description of the 
screening , properly disentangle the neutral and singly ionized Vq from the conduction band . 
According to hybrid functional calculations, all vacancies are stabilized and Vq with 2+ charge 
state is the most stable vacancy in the whole Fermi-level range inside the gap and hence Vq acts 
as a shallow donor. 

PACS numbers: 
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Ti02 has attracted many interests with apphcation as photo catalysis, water sphtting, 
solar cell and sensors^ It has been demonstrated that Ti02 can be reduced easily and reach 
high degree of nonstoichiometry which leads to high n-type conductivity^. With thermal 
excitation, oxygen vacancy (Vo) niay loose its electrons and becomes charged. Some ther- 
mogravimetric measurements reported that the electrical conductivity demonstrate various 
scaling as a function of oxygen partial pressure and temperature . Each scaling behavior 
correspond to the dominance of solely one vacancy charge state-. On the other hand , there 
are experimental findings that shows rather dominance of Ti interstitial and independent 
picture of conductivity in terms of oxygen partial pressure^. On the theoretical side, several 
investigations have been performed from embedded cluster to ab-initio total energy methods 
on the oxygen-deficient rutile Ti02 sample.-"-. Cho et al.- using local-density approximation 
(LDA) observed the entangled Vq state inside conduction band. Park et al.- in an alterna- 
tive ab-initio study using Hubbard correction for both Ti — 3d and O — 2p were only able 
to locate Vq below conduction band minimum about 1 eV . But it doesn't give the neutral 
vacancy Vq state in the gap. Moreover both above mentioned theoretical studies suffer from 
the band gap underestimation. In spite of wide research on rutile including the role of Vq in 
n-type electrical conductivity, its localized nature is still under debate. Therefore it seems 
necessary to reinvestigate the features of Vq in Ti02- Although Density Functional Theory 
has been regarded as a valuable tool for microscopic understanding of defect mechanisms, 
but the inadequate description of self-energy and underestimation of energy gap in popular 
local density (LDA) and generalized gradient approximation (GGA) leads to a widespread 
deficiencies. These limitations lead to artificial characterization of defect states and wrong 
position of thermodynamic transition levels. It is believed that overcoming gap problem 
leads to reliable determination of defect states motivated using methods such as DFT+f/, 
Hybrid functionals and Self-interaction correction to partly improve the self-interaction. 
Qualitatively, Many Body Perturbation Theory (MBPT) in the GW approach has shown 
successful treatment of quasi-particles band structure in semiconductors and is considered 
as a promising method to study defects^. However, application of GW on systems with 
partially filled d or f orbital due to wrong-predicted LDA band ordering yield unreliable 
electronic spectra and the dependence of spectra on adjustable parameters or starting point 
remain involved^. Alternatively, LDA or GGA starting point might be substituted with one 
giving closer features to respective QP band energies. As a benchmark, we choose GGA-|-?7 
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single-particle solution as starting point in GW scheme to evaluate problematic positioning 
of Vo in reduced Ti02-x- Moreover, this choice may provide insightful information about 
localized nature of the Vq states. Our results show smoothly [/-dependence of vacancy 
affinity all are grater than CBM affinity and therefore GqWq%GGA+U fails to stabilize Vq 
and Vq inside the gap. Additionally, in the realm of post-DFT methods, we employ Hybrid 
functionals PBEO and HSE in which they partially reduce self-interaction errors and improve 
gap in the description of vacancy states. The obtained results show that Vq acts as a shallow 
donor with Vq'^ being the most stable charge in the whole variation of Fermi level. 

I. COMPUTATIONAL DETAILS 

The calculations are performed using GGA+f/ based on PBE formalism for general- 
ized gradient approximation^^ as starting point for GqWq and PBEO and HSE06 exchange- 
correlation functional with a plane wave basis set and norm-conserving pseudopotentials. 
Atomic orbitals 3s, 3p, 4s, and 3d are included into the Ti valence subspace, while 3s and 
3p orbitals are considered as valence for O. Spin-polarized calculation are considered for 
systems with an odd number of electrons. The cutoff energy 70 Ryd was chosen based on 
convergence of the cohesive energy of rutile Ti02 and for Brillouin zone integration we use 
a mesh of 2 x 2 x 2. To remove the interaction of neighboring supercell we performed 72, 
96 and 108 atoms supercell calculations in which the difference between 72 and 96 atoms 
supercell defect formation energy is less than 0.02 eV and 0.04 eV for neutral and charged 
vacancies and therefore choose 72 atoms supercell. All defected supercells are fully relaxed 
until the forces are below 1 mRyd/bohr. In GGA+f/, the value of the Hubbard term U 
, applied to the Ti 3d states, is varied between and 4 eV. This range was chosen based 
on the estimate of the value of U obtained using a constrained GGA calculation, as imple- 
mented in the QUANTUM ESPRESSO codei^,and also by calculating U for an isolated Ti 
atom, divided by the theoretical optical dielectric constant eoo= 6.7. Both approaches give 
similar U values of about 1.2 eV. The GW approach is applied non-self-consistently {GqWo 
), exploiting the first order expansion of the self-energy^>^. The dynamic behavior of the 
dielectric matrix is determined by the plasmon-pole approximation as implemented in the 
SAX codei^. The defect formation energies are calculated according to 

Ef [Xg] = Etot [X,] - Etot [bulk] +^Q + q[Ep + Evbm + AV] (1) 
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where i^totf-^g] and E'totfbulk] are defected and perfect supercell total energy, respectively. 
The reference for the oxygen chemical potential is chosen to be I/2EO2 , while for Ti the 
chemical potential is referenced to the bulk Ti . The electronic chemical potential(Fermi 
energy) is referenced to the valence band maximum (VBM) corrected by the alignment 
of the electrostatic potential in perfect and defect supercell far enough from defect states. 
Moreover ,in Charged cases, the Makov-Payne correction=S have not been included due to 
large dielectric constant entering in Madelung term. 

II. GoWo@GGA+U METHODOLOGY 

It was argued that GGA+U can be viewed as an approximation to GW for localized 
d and / states^^. However , due to U correction on localized states,the hybridization of 
localized states with others has not been included i.e. the itinerant states remain at GGA 
level that in some cases isn't true . Moreover , in GGA+f/, the double counting (in the 
LDA+f/ formalism , double counting removes a part of electron-electron interaction that 
was already included in LDA Hamiltonian) term is not well defined. On the other hand 
screening is described statically while in fact it behaves dynamically and has strong energy 
dependence for particularly localized electrons . Furthermore it was shown that for open- 
shell and shallow d systems the application of GqWq on GGA doesn't yield fine QP band 
structure attributed to wrong GGA-derived band ordering^^. One solution here is to replace 
the GGA with GGA+f/ method which its eigenvalues and eigenfunctions are relatively closer 
to QP electronic structure^i^. Despite GGA+f/, GqWo calculation based on GGA+f/ has 
not shown f/-dependence explicitly while the resulted QP state would show implicit U- 
dependence^S. In the formalism of GGA+f/ based GqWq^^, the QP equation is expressed 
as: 

^ e^l' + Zi < - v., - Vdb|0f^ > (2) 

which Vdb accounts for double counting term. The GGA+f/ single particle equation is: 

eT =< - + Vgga + >= ef^+ < <p'!'\VM' > (3) 

where q can be considered as GGA energies obtained with GGA+f/ wave function. Even- 
tually, with linear expansion of self- interaction operator S(e^*) around q the QP equation 
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can be written as: 

= + Z: < <Pl'\E{e^n ' V.c\<P^' > (4) 

which clearly demonstrate advantage of having no double counting term . Regarding the 
first term, U correction affects not only the energies of the levels, but also the hybridization 
between the orbitals, if allowed by symmetry. The most pronounced modification are ex- 
pected to appear near the Fermi level, where a significant change in localization of electronic 
states can occur. The states deeper in energy are already quite localized and atomic-like, 
so that a non-zero U would mostly result to a shift in their energies. In this case, the lo- 
calized d or f occupied (unoccupied) states are pushed toward lower (higher) energies and 
get more localized. But the situation for half-filled states are quite demanding and severely 
depends on the nontrivial hybridization of included orbital as a function of U and it's char- 
acter. The second term would results in screening alteration and get tendency to increase 
energy gap with U addition. In general both terms of above equation have contributed into 
[/-dependence of the system in GW@GGA+U scheme. 

III. RESULTS AND DISCUSSION 
A. A.GW@GGA+U 

In previous section , we reviewed the application of GW@GGA+U to determine the 
affinity and ionization energy i.e. band gap of the perfect system in which the CBM is fully 
empty. In the case of defected sample which defect-induced state lies inside the gap or in 
resonance with band edges, the situation slightly changed so that the highest occupied states 
might be the defect state (in the case the defect states falls inside the gap) or CBM (in which 
the defect states is in resonance or above CB). It was shown that for Ti02 all charged Vq 
states becomes merged in CB in CCA or even GGA+f/ derived band structure. However, 
the Motivation here for the choice of GGA+f/ is that the application of U would decreases 
the outward relaxation of vacancy neighbors and enhance their hybridization leads to better 
description of the defected sample band structure closer to its real quasiparticle spectra. 

Prior to studying defected material, it is advantageous to inspect GoWo@GGA+U 
electronic structure of the perfect bulk. FigH] shows comparison between GGA+U, 
GoWo@GGA+U and XAS experimental density of states of bulk Ti02^^- Experimental 
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XAS and GGA+U was aligned at upper valence band edge with GqWo spectrum. The 
sharp peak of empty Ti — 3d in XAS spectra is situated between its analogous in GGA+U 
and GqWq . While the addition of U would enlarge the GGA+U gap and the 3d peak of Ti 
becomes closer to XAS counterpart, GqWq gives overestimated gap with 3d peak gradually 
goes away from XAS correspondent. 




energy(eV) 

FIG. 1: The Dos of Ti02 from GGA+U, GoWo@GGA+U and XAS are compared within the range 
of [/ = 0,4 . 

In both GGA+U and GqWq the most significant correction occurs for CB that are made 
of 3d-Ti via reduction of the conduction band width toward higher energies . As it is shown 
in Figj2]the [/-sensitivity is more pronounced in GqWq than GGA+U. The reason could be 
expressed in response to the question that how much the effect of U change single particle 
wave-functions and screening in dressed potential that consequently reflect the degree of 
p — d coupling . It realize that the hybridization of localized orbital with itinerant states 
that is neglected in GGA+U, play an important role when combined with GqWq . In other 
word, the differences in sensitivity of energy gap with U values could be traced back to 
the augmentation of localized and itinerant hybridization as a function of U. From the 
inspection of p-d coupling it is evident that smooth behavior of [/-sensitivity of the gap 
confirm the weakly correlated nature of the Ti02 that already was obtained about 1.2 eV in 



6 



previous section. Another useful criteria derives from U dependence of electronic structure 
is the dimensionless U/W ratio that shows the degree of the localization in Ti02 typically 
varies between 0.5-2 for the selected range of U = 1,2,3,4(W is the 3d — Ti bandwidth) . 
By comparison with the typical U/W values obtained for highly correlated material such as 
NiO which exceed from 4M , it illustrate the weak localization of d-electrons in TiOo. 




FIG. 2: The values of band gap (left) and the p-d peak distance are shown as a function of U from 
GGK+U and GoVFo@GGA+f/. 

In the case of defect involved Ti02, with the addition of U, GGA+f/ resulted Vq and 
Vq"*" energy states remain within the conduction band whereas for U = 2,3,4 , Vq state 
appear below CBM with the values 0.73, 1.12 and 1.73 eV , respectively. Figj3] shows the 
expected downward shift of fully occupied CBM in Vq and half-filled defect states in Vq 
and upward shift of unoccupied CBM in Vq and Vq~^ though the gap problem exist. Having 
earned the Vq states of GGA+f/, we are able to evaluate the influence of hubbard correction 
U compared to GW@GGA+U one. 

The GoWo@GGA+U corrected quasi-particle states are shown for all three charged states 
of Vo relative to initial GGA+f/ in FigJH Since the most GW codes does not allow a direct 
treatment of systems with unpaired electrons, the Vq defect affinity is estimated using the 
relation A(+/2+) = I (2 + /+) where A and I are affinity and ionization of given defect 
state for local charge addition. For this reason, we have performed a GqWq calculation with 
doubly ionized vacancy at the geometry of relaxed Vq . Calculated configuration coordinate 
diagram is shown in FigEl also clearly indicate almost closer configuration of Vq'^ and Vq 
and proximate same values of I(+/2+) at -R+ and R2+ ■ The affinity of the whole of the 
system was measured relative to upper valence band edge VBM while the affinity of defect 
states was measured with respect to the lower conduction band CBM. For each Vq, GqWq 
corrected states shifts up within the conduction band. 




energy{eV) 



FIG. 3: The effect of U correction on the electronic structure of perfect bulk and different charged 
oxygen vacancy is compared. For Vq with U addition the defect induced state appear in below 
CBM. 

As a result, it is observed that empty Vq'^ and half-filled Vq is more sensitive to the 
self-interaction correction than filled Vq and the shift of their quasiparticle energies is more 
remarkable than Vq. Fig. [5] shows the affinity of CBM and vacancy state of Vq and Vq'^ 
species. In GoWo@GGA+U , the affinity of CBM .i.e band gap for Vq and Vq'^ decrease 
with increasing U within 0.35 eV and 0.20 eV while this values amount to 0.10 eV and 
0.15 eV in GGA+U, respectively. For the vacancy affinity the resulted values for Vq and 
V"(^+ varies within 0.16 eV and 0.48 eV for Go Wo and 0.18 eV and 0.13 eV for GGA+f/, 
respectively. Compared to GGA+f/ , it is apparent that the U dependence of GqIVq is 
slightly more for the gap and vacancy affinity . Moreover, the affinity of the oxygen vacancy 
in both GGA+f/ and Go^^o shows similar linear ascending trend whereas the trend for band 
gap show inconsistent behavior. Inspection of the quasi-particle correction, shows that the 
upward shift is more prominent for Vq^ than Vq . This indication illustrate the different 
characters for vacancy level in which the GqWq application effectively change those states 
concerned vacancy center with no electrons(with no cost of electrostatic energy). At U=2 
both vacancies involved concrete change in quasi-particle states. This , however is an artifact 
of GGA+f/ that shows new band ordering and the vacancy level crossover some states due 
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FIG. 4: Single particle GGA+U (left) and quasi-particle energies (right) for symmetric defect state 
of various Vq are shown. The self-energy correction would shift up all defect states specially singly 
and doubly ionized vacancies. 

to change in its hybridization.. 




U(eV) U(eV) 



FIG. 5: (left) and (right) represent the comparison of affinity of the whole system and the defect 
state affinity, for Vq (top) and Vq'^ (bottom) , respectively. 

Using GW correction scheme was proposed by Rinke et af^ based on separation of for- 
mation energy into local electron addition A and lattice relaxation part A, we obtained the 
corrected formation energy and thermodynamic transition levels. At the first step, the ver- 
tical electron affinity between two charged states at fixed geometry is calculated by GqWq 
formalism . In a second step, the lattice relaxation are measured after charge addition 
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between initial and final geometries with GGA+f/.(Fig |6]) 




0.5 1 1.5 

Configuration Coordinate (q) 



FIG. 6: The simple calculated configuration coordinate diagram with GGA method. Each charged 
vacancy was computed at three configuration {Rq, i?+ andi?2+)- The diagram show closer values 
of and /_|- for singly and doubly ionized vacancies. 

Starting from reference formation energy of Vq~^, the corrected formation energies of Vq 
and Vq are then given by 

Ef{+, ej) = A(+, R+, R2+) + A{2+, R2+) + Ef{2+, = 0) (5) 

and 

Ef{0, ej) = A(+, RO, R+) + A{+, RO) + Ef{+,ef = 0) (6) 

Using calculated affinities for all defects involved we calculated formation energies as a 
function of U summarized in tableHl 

The resulted formation energy for Vq varies within 1.2-2.5 eV and within the rang of 
0.15-0.95 eV for Vq as a function of U. The increase of Vq formation energy relative to 
neutral ones implies that its concentration is now much lower than the GGA+U counterpart 
as is shown in FiglT] . Also affinity of Vq changes inasmuch as its stability as a variation 
of Fermi level (chemical potential) in the gap is vanished relative to GGA+f/ . With U 
addition, this trend will continue up along with the Vq"*" becomes the most covered stable 
charge in the energy gap and consequently the charged transition level e(2 + /O) would shift 
up though remain inside the gap. Furthermore, the charged transition level e(2 + /O) change 
from 1.70 and 1.89 eV at = 1 up to 1.82 and 2.28 eV at t/ = 4 for GGA+f/ and GqWo 
, Respectively. Hence, although GoWq@GGA+U doesn't stabilize Vq and Vq in the gap,it 
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TABLE I: All corrected formation energies by GW@GGA+U are listed in units of eV. The starting 
formation energy of V^^ has been taken from GGA+C/ values. 

A2 Ai A(+/0) A(2+/+) Efi2+) Ef{+) Ef{0) 



GGA+U GOWO GGA+U GOWO 



u= 


=0 0.05 0.01 


0.75 


2.66 


0.66 


2.16 


3.37 


3.98 


4.14 


u= 


=1 0.25 0.04 


0.79 


2.71 


0.85 


2.33 


3.81 


4.25 


4.64 


u= 


--2 0.36 0.08 


0.67 


2.43 


1.00 


2.41 


3.79 


4.53 


4.54 


u= 


=3 0.47 0.10 


0.82 


2.87 


1.22 


2.42 


4.56 


4.88 


5.48 


u= 


=4 0.55 0.15 


0.90 


2.96 


1.42 


2.37 


4.93 


5.02 


5.98 



gives Vq'^ the most dominant stable vacancy with growing U and eliminate the stability of 
singly ionized Vq for all fermi-level position within the gap. 

B. B. Hybrid functionals, HSE06 and PBEO 

We now turn to hybrid functionals that hopefully overcome the band gap problem. While 
a number of hybrid functional have been introduced^, we already focus on PBEO and 
HSE06 which typical gap has qualitative agreement with experiment. While recently pub- 
lished study with HSE06 functional gives comprehensive description of oxygen vacancies^i, 
the comparison of HSE06 and PBEO results gives efficient insight how does the short and 
long range of exchange interaction affect the defect levels. In principle ,for both PBEO and 
HSE06 ,the exchange energy part is taken as 

E^^y = aE^^ + (1 - a)E^''^ (7) 

with a=0.25 and the correlation part remains at GGA level except that in the latter the 
short-range part of the HP exchange potential is kept^. 

The reason behind choosing these functionals lies in the fact that they provide improved 
gap in many materials^. The calculations performed in this section was done with same 
parameter and accuracy as in previous section. Purthermore, our obtained total energy 
include singularity correction of the G=0 exchange potential . Por the perfect Ti02 bulk 
the energy gap are 4.15 eV and 3.10 eV for PBEO and HSE06 ,respectively shows good 
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FIG. 7: The comparison between GGA+C/ and correction scheme using GoM^o@GGA+?7 is shown. 
GqWq correction method results in metastable Yq with variation of Fermi energy. 
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FIG. 8: The DOS of hybrid functionals HSE06 and PBEO compared to GGA and GGA+[/ are 
plotted. A portion of HF exchange addition in hybrid functionals give raise to gap opening. 

agreement of HSE06 gap with experimental value 3.05 eV. The comparison of different 
functional gap are shown in FiglHl from inspection of above density of states it seems that 
the addition of a portion of HF exchange would considerably increase the gap relative to 
GGA. The comparison of PBEO and HSE06 also reflect the fact that the energy gap is 
correctly reproduce by the short range nature of exact exchange. The VBM state lowering 
and GEM state lifting amount to 0.96 eV and 1.19 eV for PBEO and 0.63 eV and 0.77 eV 
for HSE06. This trends is attributed to the HF admixed portion which reduce the self- 
energy of VBM and also raise the GBM. Using above methods we calculated the energetic, 
electronic and structural properties of each vacancy type. Since the vacancy induced state 
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TABLE II: The magnitude of vacancy neighbors outward relaxation in A for PBEO and IISE06 are 
hsted. The relative comparison shows lower relaxation of fully and half occupied vacancy levels(due 
to bounded electron to defect center) in hybrid functionals than GGA. 



Ti-0 VS 



o 



y2+ 



GGA 1.128 1.151 1.127 1.153 1.126 1.157 
PBEO 1.014 1.120 1.042 1.076 1.130 1.163 
HSE06 1.000 1.121 1.009 1.124 1.149 1.175 




energy(eV) 



energy(eV) 



FIG. 9: Both HSE06 and PBEO are successful in giving the neutral and singly ionized vacancies 
below CBM . HSE06 band gap is coincide with experimental gap and also the defect state is in 
closer agreement with experimental values and Re^. 



are made from the three Ti atoms surrounding the vacancy , their relaxation magnitude is 
determinant to the positioning of vacancy states in the band structure. Table. HTVepresent 
the distances of vacancy neighbors for GGA, PBEO and HSE06. 

It is obvious that in PBEO and HSE06 the outward relaxation for Vq and Vq that have 
bounded electrons to vacancy center is relatively small compared to GGA while bigger 
outward relaxation is observed for Vq^ than GGA analogous due to having no bounded 
electron. Having earned the amount of vacancy relaxation, we are able to interpret the 
obtained states in the gap as shown in FigjOl 

Upon introducing oxygen vacancy , its neighbors relax in order to reinforce their bonding 
with the rest of the lattice and results in less overlap between three Ti atoms and shift 
the vacancy state up. On the other hand if the vacancy center has bounded electrons it 
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FIG. 10: The diagram of formation energy with Fermi energy variation is depicted for PBEO and 
HSE06 hybrid functional for 0-rich and Ti-rich conditions. Both functionals give the V^"*" the 
stable vacancy with Fermi energy change and transition level e(2 + /O) is approximately within ~ 
1 eV above conduction band minimum. 

resist against outward relaxation since it gives unstable energetic picture. Hence balance of 
these two competitive factors determine the position of vacancy states. From the table [Til 
it evidences that HSE06 and PBEO show less outward relaxation than GGA for Vq and Vq 
while the situation is reversed for Vq'^ which give raise to the Vq and Vq levels 1.12 and 
1.25 eV for PBEO and 0.80 and 1.40 eV for HSE06 below CBM . Comparing Vq states for 
neutral and singly ionized charged vacancy in PBEO and HSE06 indicate that energy shift 
of Vq state is more than Vq due to the presence of short part of exact exchange. Finally 
we address the extrapolated formation energy at 0-rich and Ti-rich limit as shown in Fig. 
[TUl The formation energy of Vq'^ is fairly low compared to Vq and Vq in terms of Fermi 
energy variation in the gap . The reason lies in the fact that these hybrid functionals lower 
the VBM level considerably lead to lowest formation energies(with no electron bounded 
to vacancy) and therefore the charge transition level lies above the gap . Our results for 
formation energy and e(2 + /O) is closed to HSE results of Ref.^^ . This implies that both Vq 
and Vq have always higher formation energy than Vq'^ make it only possible stable charged 
state of oxygen vacancy. This behavior is exactly same in both 0-rich and Ti-rich regime 
except that the formation energies in Ti-rich are much lower than 0-rich. Accordingly, 
this observation is mainly consistent with experimental measurement suggest conductivity 
reduction with increase of oxygen partial pressure^i. 
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IV. SUMMARY 



In conclusion, we studied the problem of oxygen vacancies in viewpoint of gap correction 
method. With the apphcation combined GW correction and GGA+U starting point referred 
as GW@GGA+U, we found that the quasi particle energies shows different [/-dependence 
traced back to their character and occupancy.. Additionally, we found that application of 
GqWq depends on the starting point . It also fails to predict the vacancy state inside the 
gap contrary to the experimental findings. Alternatively , we employed hybrid functional 
PBEO and HSE06 in which they are able to truly predict the neutral and singly ionized 
vacancies in the energy gap . However they gives raise to stable energetic picture of doubly 
ionized vacancy as a function of Fermi energy variation within the gap and transition level 
e(2 + /O) about 0.6 eV and 0.8 eV for PBEO and HSE06 above CBM reflect shallow n-type 
conductivity features of Vq- 
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